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We examine the effect of short unstable periodic orbits 
on wavefunction statistics in a classically chaotic system, and 
find that the tail of the wavefunction intensity distribution in 
phase space is dominated by scarring associated with the least 
unstable periodic orbits. In an ensemble average over systems 
with classical orbits of different instabilities, a power law tail 
is found, in sharp contrast to the exponential prediction of 
random matrix theory. The calculations are compared with 
numerical data and quantitative agreement is obtained. 

Quantum eigenstates of classically chaotic systems 
generically exhibit a phenomenon known as scarring, the 
enhancement of intensity along short unstable periodic 
orbits for some fraction of the wavefunctions. Scarring is 
a fascinating example of the influence of identifiable clas- 
sical structures on stationary quantum properties and on 
long-time quantum transport in a classically ergodic sys- 
tem. The occurrence of scars is in some sense paradoxi- 
cal, because classically, all such short-time information is 
destroyed at long times, and a classical probability dis- 
tribution after being evolved for a sufficiently long time 
retains no memory of its initial state. Scarring is one of 
the most dramatic examples of a departure of quantum 
chaotic systems from the predictions of random matrix 
theory (RMT), according to which wavefunctions must 
be evenly distributed over phase space, up to quantum 
fluctuations. Scarring has now been observed experimen- 
tally in a variety of systems, including microwave cavities 
jl],|]], tunnel junctions [3], and the hydrogen atom in a 
uniform magnetic field \4 |J . 

Examples of scarring were observed numerically in |jj , 
and a theory based on the semiclassical evolution of 
Husimi states near a periodic orbit was provided. Later 
work by Bogomolny M and Berry || involved calcula- 
tions in coordinate space and Wigner phase space, re- 
spectively. All these works were based on the linearized 
dynamics around the unstable periodic orbit, and were 
thus, by construction, theories of the short-time behav- 
ior only. Yet to get a true understanding of the proper- 
ties of individual eigenstates it is essential to understand 
the long-time quantum dynamics, including returns of 
amplitude to the original periodic orbit after undergoing 
excursions into other areas of phase space. In a recent pa- 
per H , a formalism was developed for dealing with these 
nonlinear contributions to scarring, providing quantita- 
tive agreement of the theory with numerical results. This 



work used a measure of scarring based on Husimi in- 
tensities. (Recently Fishman, Agam, and others have 
provided interesting new perspectives on the problem on 
scarring, and have offered a measure of scarring related 
to, but somewhat different from ours [fTof . A number of 
other authors have also made significant contributions in 
this area; we cannot list them all but a few recent refer- 
ences are provided in pd[|.) In this paper, we will apply 
a result previously obtained in j^] to compute the effect 
of scarring on the wavefunction intensity distribution, a 
function which has been investigated previously for dif- 
fusive systems [fl2"f . 

We first give a general idea of the formalism and state 
the key result of |9|] which is relevant to the present work. 
A gaussian wavepacket | W ) is initially centered close to 
a classical periodic orbit and allowed to evolve. If the 
instability exponent A of the orbit is small, the func- 
tion A(t) = (^>\^>(t)) will (under the linearized dynam- 
ics) exhibit recurrences on a time scale of one period 
Tp, and the amplitude of these recurrences decays on 
a scale Tp/X. This decay leads to the formation of en- 
velopes in the local density of states S(E), which is the 
fourier transform of the autocorrelation function. These 
envelopes in the energy spectrum have spacing h/Tp, 
width scaling as hX/Tp for small A, and height scaling 
as A . Long-time (nonlinear) recurrences lead to fluc- 
tuations multiplying this envelope in the energy domain. 
Eventually, by the Heisenberg time Th = S-/A, where A 
is the mean level spacing, individual peaks are resolved 
in the spectrum S(E), the heights of the peaks being the 
intensities of the corresponding eigenstates at the test 
state |*), S(E) = \(n\V)\ 2 5(E - E n ). 

The simplest picture of the non-linear recurrences as- 
sumes that the periods, actions, and homoclinic points 
corresponding to the long-time excursions are random 
and uncorrelated, up to the constraint of unitarity which 
produces delta-function spectral peaks by the Heisenberg 
time. This, when combined with the known short-time 
dynamics in the linear regime, can be shown to produce a 
chi-squared distribution of spectral intensities multiply- 
ing the original linear envelope]^] 
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1 Strong, isolated nonlinear recurrences which cannot be 
treated statistically, as well as correlations between long-time 
excursions modify this simple picture. However, we will find 
that for the generalized baker's maps (a paradigmatic exam- 
ple of hard chaos), the assumptions of randomness beyond the 
linear decay time lead to results which are in quite reasonable 
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A key result of the calculations in Ref. |^ is the follow- 
ing: individual spectral lines (overlap intensities) in the 
local density of states obey the usual chi-squared (Porter- 
Thomas) fluctuations, but these are modulated by the 
fourier transform of the linearized short-time autocorrela- 
tion function. The latter can be computed analytically in 
terms of the instability A of the classical orbit in question 
(see Eq. |] below) . More explicitly, in the case of complex 
eigenstates, the chi-squared distribution has two degrees 
of freedom, and in the absence of scarring the probability 
of having a spectral line height greater than x is given by 
P{x) — exp(— x). Here x is normalized to have a mean 
value of unity, i.e. x n = N\(n\^)\ 2 , where TV is the total 
number of states. Now in the presence of scarring this is 
modified to 



P(E,x)=exp(-x/S iia (E)) 



(1) 



for a eigenstate with energy E. Here Su n is the spectral 
envelope given by the fourier transform of the linearized 
dynamics (Eq. ^ for the case treated in the present pa- 
per). 

In the following discussion we will for simplicity con- 
sider the case of scarring by a fixed point of a discrete- 
time map. (The generalization of the results to orbits of 
period greater than one and to continuous time systems 
is straightforward ||.) So let us consider without loss 
of generality a hyperbolic fixed point at the origin of a 
compact phase space, with exponent A, and stable and 
unstable manifolds oriented along the p and q axes, re- 
spectively. The equations of motion near the fixed point 
are then given by 



q' = gexp(At) 
p' = pcxp(— At) 



(2) 



In the presence of shearing, e.g. dq' /dp ^ 0, or for non- 
orthogonal stable and unstable manifolds, or in the case 
where the manifolds are not oriented along the p and q 
axes, a canonical transformation would first need to be 
performed to get the equations of motion into the form 
above. Now we define our test state to be a gaussian 
wavepacket centered at (qo,Po), with horizontal width a 
and vertical width o~ p = %/a. In coordinate representa- 
tion this is given by 



*(?) = (4iral) 1/4 exp[-(q - q ) 2 /2<r 2 + ip (q 



qo)/h] . 
(3) 



In situations described above where the local equations 
of motion do not have the form of Eq. || in the natural 
coordinates, an optimal test state would have a complex 
width a in those coordinates, as can be seen by perform- 
ing a canonical transformation of the gaussian of Eq. H 



agreement with numerical data. 



We now allow the wavepacket to evolve under the lin- 
earized dynamics, stretching each time step by a factor 
e A in the q— direction and shrinking by the same factor in 
the p— direction. The linearized quantum autocorrelation 
function after time t is given by 



A U n(q,P,cr,X,t) = 



exp i8ot 



x exp 



cosh \t 
cosh Xt — 1 



2 cosh Xt \ a 



iqp 



tanh Xt 



(4) 



Here 9q is the phase associated with one iteration of the 
periodic orbit, given by the classical action in units of 
H, plus any Maslov indices associated with caustics in 
the classical dynamics. 8q determines the location of the 
peak in the spectral envelope Su n (E), defined to be the 
fourier transform of An n (t). Since we will be interested 
in performing an energy average, we freely set 9q = 0. 

Now the expression in Eq. |] can be inserted into Eq. |l| 
to obtain the distribution of wavefunction intensities at 
a given energy. Because we wish here to consider all 
eigenstates, independent of energy, we then perform an 
energy averaging, remembering that for a map the quasi- 
energy is defined to lie between and 2ir only. We also 
notice that the tail of the intensity distribution will be 
dominated by the peak of the spectral envelope at E = 0, 
and we therefore can use a saddle point approximation, 
obtaining 



P(q,p,a,X,x) = — J dEP(q,p,a,X,E,x) 



1 exp(-x/ S\j n (q,p, a, A, E = 0)) 

-x d 2 S lin 
5 lin (0)^ dE* 



2tt 



(5) 



where the expression obtained is an asymptotic form 
valid for large x. For small A, the sum over time steps can 
be replaced by an integral, and we have at qo = po = 



S lin (E) = / dt 



iEl 



V cosh Xt 



(6) 



analysis, Sii n (0) = Q/X and 
where Q and W are numerical con- 



Now by dimensional 

&(0) = -W7A 3 
stants. We thus obtain the first result of this paper, the 
tail of the intensity distribution for a wavepacket centered 
on a periodic orbit, 



P(q = 0,p = 0, a, X, x) 



1 



Q 



2tt VW 



X{xX) 



-l/2 e -x\/Q 



(7) 



Notice that the exponential tail has been effectively 
stretched by a factor of Q/X, corresponding to the in- 
creased height of the peak of the linear envelope at small 
A. There is also a linear suppression factor of A, corre- 
sponding to the width of the peak in Su n (E), and indicat- 
ing that only a fraction scaling as A of all the eigenstates 
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are effectively scarred. The remainder of the eigenstates 
are typically "antiscarred" , having on average a lower 
intensity at the periodic orbit than would be expected 
based on RMT. This distribution will have a nontrivial 
inverse participation ratio (IPR), scaling as the inverse 
of the width of the linear envelope, i.e. as 1/A. 
The region of validity of Eq. fj] is 



l«A" 1 «i«JV. 



(8) 



The first inequality ensures that many iterations of the 
periodic orbit contribute (so the sum over iterations can 
be replaced by an integral) and the scarring is strong. In 
fact, however, because of the large value of the numerical 
constant Q, the formula works well even for exponents 
as small as log 2, as will be seen in the numerical study 
below. The second inequality says that we are in the tail 
of the distribution and the events are all coming from the 
peak of the linear envelope. The third inequality is a uni- 
tarity constraint - obviously our assumption of random 
fluctuations breaks down for intensities of order unity, 
when the entire wavefunction would be concentrated in 
a phase space area of order %. 

Now we go on to perform a similar analysis integrat- 
ing over the phase space variables qo and po- As be- 
fore, we take the exponential exp (— x/ Su n (qo,po, . . .)) 
and expand to second order in qo,po around the max- 
imum go — po = 0. Then upon integration by stationary 
phase we pick up a determinant factor of 



27rS lin (0) 2 



1 



d 2 5l in d 2 Sl in 

dql dpi 



(9) 



Here we have taken the classical phase space volume to 
be unity for simplicity. Now for small A, 



<9 2 <Slin ,„n _ 2 — Z 



A 



(10) 



where Z is yet another numerical constant, and similarly 
for po, with <7 replaced by a p =fi/a. So the total factor re- 
sulting from the phase space integration is ^- , again 
independent of a. Combining this with the expression in 
Eq. [7] above, we obtain the second result, for the distri- 
bution of overlap intensities after energy and phase space 
averaging, 



V Q 3 



2 ZVW 



hX(x\)-^ 2 e^(-xX/Q). (11) 



Here we have picked up a factor of U from the factor of 
(7 in Eq. [l^ and the corresponding factor of <j p — h/a 
associated with the falloff in S\ ln in the momentum di- 
rection. This indicates that the tail is coming entirely 
from the region near the periodic orbit, specifically from 
wavepackets that have large classical probability density 
right on the orbit. (Thus, a measure like the IPR for a 
generically placed wavepacket will not see the effect of 



scarring by an individual periodic orbit, when the semi- 
classical limit has been taken.) The result in Eq. [ll] is 
valid in the regime 



max(log N, X l ) < x < N . 



(12) 



Here log N is the value of x near which the RMT expo- 
nential decay law reaches values of order h — l/2irN. In 
this region, a crossover occurs between the head of the 
distribution, which is dominated by non-scarred region of 
phase space and approaches the Porter-Thomas (RMT) 
prediction, and the tail, dominated by scarring, given by 
the expression above. The expression Eq. [ll] holds also 
for an ensemble of systems, all having an orbit with insta- 
bility A. In principle, we should of course do a sum over 
all periodic orbits, however the tail will clearly always be 
dominated by the orbit with smallest A. 

Finally, we now consider an ensemble of systems where 
the value of the smallest exponent varies from system to 
system, with distribution V(\) = CA Q for small A. Then 
using Eq. [ll] and integrating over A we obtain 



JlJ^Q^/2 T{a + l/2)hx-^+ a y (13) 

V 2 zVw 



P(x) = C 



Note that this is an uncontrolled approximation because 
we have integrated over A after having assumed xX was 
large. However, if we had included higher-order correc- 
tions in (xX)^ 1 in Eq. [ll], the scaling of P{x) would re- 
main unchanged, i.e. 



P{x) = Cf(a)hx' ( - 2+a) , 



(14) 



with the dimensionlcss function /(a) somewhat differ- 
ent from that given in Eq. |l3|. An important point is 
that the tail displays power-law behavior in the intensity 
x, a strong deviation from the exponential prediction of 
RMT. As with Eq. [ll], this asymptotic form is valid for 
values of x large compared to log N and small compared 
to N. For small x we again expect a crossover to the 
Porter-Thomas form. For large x we expect a downward 
correction away from the ;£-( 2 + Q ) form, with a break- 
down of the approximation occurring at some fraction of 
N, depending on a. 

Now, we proceed to test numerically these predictions 
of the nonlinear scarring theory. The system we use for 
this purpose is the generalized three-strip baker's map, 
described in some detail in Q . This system is similar to 
the original baker's map, except that the two strips are 
replaced by three, with widths generally unequal, but 
normalized to X)i=o w * = There is a fixed point of 
the classical dynamics associated with the middle strip, 
and the instability associated with this orbit is given by 
A = I log iui I . So we choose wi — 1/2, set N — 200, 
and find numerically the wavefunction intensity distribu- 
tion at the fixed point after ensemble averaging over the 
widths wo, 2- (The predictions are expected to hold for 
individual systems as well, at sufficiently large values of 
N. However, for the matrices which we can efficiently 
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diagonalize, N is not large enough to obtain good statis- 
tics in the tail without ensemble averaging.) A circular 
wavepacket of width a = \f% is used. The results are 
compared in Fig. [I] (upper thick curve) with the pre- 
diction of Eq. [?], plotted as a dashed curve. The RMT 
prediction is shown as a dotted line. Note that the theo- 
retical prediction of Eq. [7] has no free parameters and is 
^-independent, depending only on the exponent A of the 
periodic orbit in question. 

Next, we perform a phase space average for the systems 
described above, collecting statistics for wavepackets uni- 
formly distributed over the entire phase space. The re- 
sulting statistics are also plotted in Fig. [I] (lower thick 
curve), where the theoretical prediction for the tail, given 
by Eq. [ll], is shown as a solid curve. Again, the Porter- 
Thomas distribution appears as a dotted line. We see 
a crossover between the two regimes at a value of x of 
order log N. 

Finally, we want to construct an ensemble which will 
contain systems with orbits of different instability ex- 
ponents A. For this purpose, we take a uniform dis- 
tribution of strip widths Wq and W2, each in the range 
[0, 1/4]. The fixed point in the middle strip, with expo- 
nent A = I log(l — wo — W2)\, is always the least unsta- 
ble periodic orbit. This ensemble has power a = 1 and 
C = 4 2 = 16 in the notation of Eq. |l3|. Averaging over 
100 systems, we obtain the statistics plotted in Fig. ||. 
The power-law prediction given by Eq. [l3] is plotted as 
a solid line on the log-log scale, with the RMT predic- 
tion as a dotted curve. Once again, we see a crossover 
between the two regimes for x of order logiV ss 5.3. We 
also see a gradual breakdown of the approximation as x 
approaches values comparable to N — 200. Note that 
the quantitative agreement is in spite of the fact that an 
uncontrolled stationary phase approximation was used 
in obtaining the overall constant in front of Eq. |l3|, as 
explained above. The important thing to notice here is 
the power-law behavior of the tail, in agreement with the 
theory, and the dramatic deviation from the predictions 
of RMT. By x — 100, where the approximation x <C N 
is clearly beginning to break down, the measured proba- 
bility is still within a factor of 4 of our prediction and is 
enhanced by 10 37 over the Porter-Thomas value. 

We have also checked the linear ?i-dependance of the 
phase-space averaged results Eqs. [ll],|l3| by repeating the 
preceding numerical analysis with larger matrices [N — 
500, 1000). In addition, we have constructed an a = 
ensemble by imposing the restriction wo = W2 and have 
observed a x~ 2 power-law behavior in accordance with 
Eq. H 



useful discussions. 
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FIG. 1. Cumulative wavefunction intensity distribution (a) 
as measured by a test state centered on a periodic orbit with 
instability A = log 2, plotted as the upper thick curve with 
scarring theory prediction given by dashed curve, and (b) 
averaged over the entire phase space of size 200/i, plotted 
as lower thick curve with theory given by solid curve. The 
dotted line is the Porter-Thomas law. 

FIG. 2. Cumulative wavefunction intensity distribution af- 
ter ensemble averaging over systems with classical orbits of 
different instability exponents. Here again N = 200, and the 
dotted curve is the RMT prediction. 
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